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Depth of investigation in electromagnetic sounding methods 



Brian R. Spies* 



ABSTRACT 

The time or frequency at which the electromagnetic 
(EM) response of a buried inhomogeneity can first be 
measured is determined by its depth of burial and the 
average conductivity of the overlying section; it is 
relatively independent of the type of source or receiver 
and their separation. The ability to make measure- 
ments at this time or frequency, however, depends on 
the sensitivity and accuracy of the instrumentation, 
the signal strength, and the ambient noise level. These 
factors affect different EM sounding systems in sur- 
prisingly different ways. 

For the magnetotelluric (MT) method, it is possible 
to detect a buried half-space under about 1.5 skin 
depths of overburden. The maximum depth of inves- 
tigation is virtually unbounded because of high signal 
strengths at low frequencies. Transient electromag- 
netic (TEM) soundings, on the other hand, have a 
limited depth of penetration, but are less affected by 
static shift errors. For TEM, a buried inhomogeneity 
can be detected under about one diffusion depth of 
overburden. For conventional near- zone sounding in 
which induced voltage is measured (impulse re- 



sponse), the depth of investigation is proportional to 
the 1/5 power of the source moment and ground 
resistivity. By contrast, if the receiver is a magnetom- 
eter (step response system), the depth of investigation 
is proportional to the 1/3 power of source moment and 
is no longer a function of resistivity. Magnetic-field 
measurements may, therefore, be superior for explo- 
ration in conductive areas such as sedimentary basins. 
Far-zone, or long-offset, TEM soundings are tradition- 
ally used for deep exploration. The depth of investi- 
gation for a voltage receiver is proportional to the 1/4 
power of source moment and resistivity and is in- 
versely proportional to the source-receiver separation. 
Magnetic-field measurements are difficult to make at 
long offsets because instrumental accuracy limits the 
measurement of the very slow decay of the magnetic 
field. 

Frequency-domain co.itrolled-source systems are 
ideally suited for sounding at the very shallow depths 
needed for engineering, archaeological, and ground- 
water applications because of the relative ease of 
extending the measurements to arbitrarily high fre- 
quencies, and also because geometric soundings can 
be made at low induction numbers. 
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INTRODUCTION 

Estimation of the depth of investigation is crucial for both 
the design and interpretation of electromagnetic (EM) 
sounding experiments. Most geophysicists know that the 
basic parameters which control the depth of investigation 
include transmitter moment, noise levels, and earth resistiv- 
ity; but the specific relationships of these parameters to the 
depth of investigation achieved are sometimes poorly under- 
stood. Survey design based on heuristic considerations or 
intuition can either fail to meet the geologic objectives or 
result in unnecessary expense. The purpose of this paper is 



to derive simple rules-of-thumb that can be used by the 
exploration geophysicist to calculate the depth of investiga- 
tion for EM depth sounding methods commonly used. 

EM sounding is widely used in petroleum, mineral, 
groundwater, and geothermal exploration and includes con- 
trolled-source techniques (both frequency- and time-domain) 
and natural-source techniques such as magnetotellurics. An 
overview of applications, methodology, and case histories 
for controlled-source techniques is given by Spies and 
Frischknecht (1989), and for magnetotellurics by Vozoff 
(1989). The techniques used most widely for geophysical 
sounding are the magnetotelluric (MT) and time-domain or 
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^sient electromagnetic (TEM) methods, which together 
comprise approximately 76 percent of EM surveys (by 
dollars spent) in the Western world (Montgomery, 1987). 
tV j agn etotellurics is used mainly for structural mapping of 
basins, especially those covered with resistive volcanic or 
overthrust material, whereas TEM is used for shallower 
sounding and stratigraphic applications. Frequency-domain 
electromagnetic (FEM) soundings are used mainly for shal- 
low engineering, archeological, environmental, and ground- 
water studies. 

It is sometimes difficult to choose the optimal EM sound- 
ing method, particularly since advances in instrumentation 
have widened the depth range beyond that which was 
possible in the past. Having chosen the method, the design 
of the survey must be such that the desired depth of 
investigation can be achieved. In some geologic environ- 
ments there may be advantages in using a combination of 
techniques. For example, corrections for MT static shifts 
based on TEM measurements (Andrieux and Wightman, 
1984; Sternberg et al., 1988) depend on an adequate overlap 
in the depths of investigation of the two techniques. 

In this paper, I relate the depth of investigation to the 
physics of EM diffusion, instrument sensitivity, and noise 
levels for the MT and TEM methods. Rules-of-thumb are 
derived which can be used by the exploration geophysicist 
both to design a field survey and to choose the optimum field 
parameters so that the EM sounding covers the desired 
range of depths. 

THEORETICAL BACKGROUND 

To understand which parameters affect the depth of pen- 
etration, it is helpful to consider the diffusion and attenuation 
of EM fields within the earth. For a uniform earth and 
spatially one-dimensional (1-D) (plane-wave) fields the or- 
thogonal electric- and magnetic-field components in the earth 
are (Nabighian and Macnae 1989) 



and 
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where E Xq is the field at the surface, and 
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Is the frequency-domain (FD) skin depth. It is assumed that 
conductivity a and magnetic permeability jx 0 are frequency- 
independent and that displacement currents can be ignored. 

s can be seen from these equations, the depth z = 5 is the 
||epth at which the EM field is reduced to l/e, or 37 percent, 
0 rts value at the surface, while the phase is rotated by one 
radian. 

T he time-domain analogs to these equations are given by 
funcr 196 ^ Nabighian and Macnae (1989). For a step- 



ton excitation of magnitude h 0 established at time / = 0, 



the transient fields 




and 



e x (z, t) = 



h y {z, t) = h 0 erfc -% / 
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(3a) 



(3b) 



where erfc is the complementary error function. For a fixed 
time f, the variation of e x (z 9 t) and h y (z, t) with depth is 
similar to the normal density and nor mal dis tribution func- 
tion with zero mean and variance V2//o> 0 . By keeping z 
fixed in equation (3a) and equating to zero the derivative of 
e x (z, t) with respect to time, the maximum transient electric 
field at any time t can be shown to be located at a depth of 




(4) 



are 



where the quantity 5 TD is the time -domain diffusion depth. 
The similarity between the expressions for diffusion depth 
[equation (4)] and skin depth [equation (2)] is remarkable; 
the depths are proportional to 1/Vco in the frequency do- 
main, and Van the time domain. As a function of time, the 
time-domain diffusion depth travels downward with a veloc- 
ity given by 



The subsurface pattern of the diffusion currents is an 
expanding, attenuating image of the source geometry, and is 
most easily visualized in the time domain. "Snapshot" 
images of the electric-field strength in a homogeneous earth 
are given by Nabighian and Macnae (1989) for a plane-wave 
source, by Oristaglio (1982) for a line source, by Lewis and 
Lee (1978), Hoversten and Morrison (1982), and Newman et 
al. (1989) for a loop source, and by Gunderson et al. (1986) 
for a grounded wire and infinite line source. 

The terms "skin depth" and "depth of penetration" are 
often assumed to be synonymous in elementary textbooks, 
although in practice it is necessary to distinguish between 
these two terms. Theoretically the EM fields at any fre- 
quency or time are present at all depths, albeit at vanishingly 
small values. The practical depth of investigation from the 
geophysical standpoint depends on such factors as the 
sensitivity and accuracy of the instrumentation, the com- 
plexity of the geologic section, and the ambient or inherent 
noise levels. In an ideal geologic environment, the practical 
depth of investigation may be several skin depths, whereas 
in geologically complex or noisy areas, and for certain EM 
sounding systems, it may be much less than one skin depth. 

Most approximate 1-D interpretation schemes (e.g., Bos- 
tick, 1977) are based on a direct transformation of the data to 
a conductivity-depth section. The ' 'effective depth of pene- 
tration*' defined in these schemes is closely related to skin 
depth or diffusion depth (see Appendix A) and is controlled 
by the conductivity of the ground and the frequency (or time) 
of measurement which, together, determine the depth to 
which the induced currents have diffused. 

The approach used in this paper is to find the highest 
frequency (or earliest time) at which a homogeneous half- 
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space buried beneath a uniform layer has a measurable effect 
on the observed response. It is argued that the expressions 
so derived are also valid for detection of 2-D and 3-D bodies 
at depth. 

DETECTABILITY OF A BURIED HALF-SPACE 

Magnetotelluric sounding 

The factors affecting the depth of investigation can most 
easily be studied by considering the suite of two-layer 
Cagniard apparent resistivity and phase curves shown in 
Figure 1. The basic features of these curves are well-known: 
the response at high frequencies is dominated by the upper 
layer, while the underlying half-space dominates at low 
frequencies. 

The "depth of investigation" is defined here as the max- 
imum depth at which a buried half-space can be detected by 
a measurement system at a particular frequency. This depth 
is dependent on the conductivity contrast and the accuracy 
of the recording and processing system. Consider a system 
having a measurement accuracy of ±3 dB in the magnitude 
of EIH, which corresponds to an accuracy of ±6 dB in 
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Fig. 1. MT Cagniard apparent resistivity and phase curves 
for a two-layer earth. The horizontal axis is the upper layer 
thickness expressed in skin depths. Dashed vertical lines 
superimposed on the figure correspond to values of 1, 1.5, 
and 2 skin depths. 



apparent resistivity. The response of a two-layer earth wj^j 
a conductivity contrast of 2 lies inside the error bounds at 
frequencies (shaded region in Figure 1). However, the re : 
sponse of the same model having a conductivity contrast of 
5 is outside the error bounds when the overburden is l 
than 0.3 skin depths thick. 

On the other hand, if the measurement accuracy we^ 
increased so that apparent resistivity could be measured to 
±0.1 dB (i.e., approximately 1 percent), then the layering 
would be observable under two skin depths of overburden 
Most interpreters, however, would not be confident i n 
basing their interpretation on a small phase shift and the start 
of a'small undershoot or overshoot in the apparent resistivity 
data (although their automated computer inversion programs 
may have no such qualms). An interpreter usually assigns a 
"comfort" or 1 'confidence' ' factor based on typical, or 
best-case, data quality and system accuracy, and expert 
ence. I assume here that apparent resistivity can be mea- 
sured with an accuracy of ±0.8 dB (approximately iq 
percent), which would enable the detection of a buried 
half-space with a conductivity contrast of 10 under 1.5 skin 
depths of overburden. Examination of other curves based on 
alternative definitions of apparent resistivity (e.g., Spies and 
Eggers, 1986) supports this observation. Furthermore; 
Szarka and Fischer (1989), by analysis of the real and 
imaginary parts of the impedance, show that the final devi- 
ation of the phase curve from 45 degrees at dlh = 1.58 
physically implies the final separation of the centers of the 
in-phase and out-of-phase current systems at depth d. Thus, 
a reasonable estimate for the depth of investigation for MTis 
taken to be approximately equal to 1.5 skin depths, or 



z* 1.5 5 FD = 1.5 



aj(JL 0 co 



= 750 




One method of estimating the contribution to the EWj 
response from various depths is to calculate the Frecheil 
derivatives from the appropriate integral equation (e.g.. 
Parker, 1977; Oldenburg, 1979; Chave, 1984; Edwards. 
1988). The Frechet derivatives give the sensitivity of thij 
measurements to small local changes in- conductivity, anij 
the shape of the Frechet kernel may be analyzed to obtaioj 
approximate penetration depths and the resolution interv; 
associated with any model and depth interval. By analysis oi 
Frechet derivatives for a plane-wave source, it is shown iij 
Appendix B that 95 percent of the MT response from af 
homogeneous earth originates from depths of less than 1 M 
skin depths. | 

.1 

Maximum depth of investigation. — There is no difficulty iflg 
extending an MT sounding to arbitrarily great depths, sincc| 
the signal strength increases approximately in proportion tc| 
\lf. The maximum depth of investigation is not limited bf| 
system sensitivity or noise levels but depends only upon thfg 
measurement frequency and earth conductivity as described 
by equation (5). There are major differences in the analogoi*| 
expressions for TEM sounding, as is shown in a lat*| ; 
section. I 

Upper limit of investigation.— The ability to resolv| 
shallow parts of the geoeiectric section is limited by *| 
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highest frequency of measurement (usually <100 Hz, but up 
to 15 kHz in some MT systems). From apparent resistivity or 
amplitude data, only the "average" conductivity cr av of the 
isoelectric section lying above the penetration depth z min 
calculated at this frequency can be determined. z min can be 
considered the minimum resolvable depth obtainable in an 
sounding and is calculated by substituting a av for in 
equation (5). The depth at which conductivity variations can 
be investigated is a function of both resistivity and frequency 
and is shown in diagrammatic form in Figure 2. 

It should be mentioned that if phase (or the derivative of 
amplitude) is available, then additional constraints can be 
placed on the conductivity structure above z min . For exam- 
ple, Schmucker (1987) describes two extremal models with 
large resistivity contrasts where the top layer acts either as a 
nonconductor (A-model) or as a thin-sheet conductor 
(T-model). For the /i-model, the phase of the impedance <(>(<*>) 
lies between 45 and 90 degrees, and for the limiting case of a 
perfectly resistive top layer of thickness h overlying a basal 
half-space, the MT response is given by 2(o))//(o = h + 
(1 - i)8 2 (w)/2, where 8 2 (w) is skin depth in the basal 
half-space. For the T-model, the phase lies between 0 and 45 
degrees, and the limiting model is that of a thin surface layer 
of conductance S overlying a half-space, for which the MT 
response is /WZ(o>) = /tu|x 0 5 + (1 + i)/8 2 (o>). If the phase is 
45 degrees, the response reverts to that of a uniform half- 
space, with resistivity equal to the Cagniard apparent resis- 
tivity. 

Depth of investigation in a multilayered earth.— The depth 
of investigation given by equation (5) can be extended to a 
multilayered earth by replacing a 1 by the "effective" or 
"average" conductivity of the ground overlying the depth of 
investigation. The efficacy of using the effective or average 
conductivity is illustrated by numerous studies of the equiv- 
alence and nonuniqueness of multilayered earth models that 
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have appeared in the inversion literature in the last two 
decades. In the case of imprecise observations (i.e., all those 
obtained in practice), individual inversion algorithms can 
yield very different solutions, and the data can be fit with a 
large variety of models. These range from the delta function 
D + models of Parker (1980), through simple layered models 
familiar to most geophysicists (Jupp and Vozoff, 1975; 
-Nabetani and Rankin, 1969; Glenn et al., 1973; Larsen,' 
1981), to smoothly varying models (Glasko et al., 1972; 
Weidelt, 1972; Oldenburg, 1979; Oldenburg et al., 1984; 
Constable et al., 1987; Srnka and Crutchfield, 1987), Exam- 
ples of vastly different models fitting an observed data set are 
given by Jones and Hutton (1979), Parker and Whaler (1981), 
Oldenburg et al. (1984), and Constable et al. (1987). 

Much more stable inversions are achieved if the solution is 
sought in terms of integrated or cumulative conductance 
down to depth z, 



5(z)= v(z)dz. 



(6) 



This result comes from the classic study of thin sheets by 
Price (1949) and has been used by numerous other workers, 
e.g., Schmucker (1970), Fullagar and Oldenburg (1984), and 
Macnae and Lamontagne (1987). Oldenburg (1983) and Wei- 
delt (1985) discuss the issue of extremal bounds on inte- 
grated conductance estimates. 

To estimate the depth of investigation in a multilayered 
earth section, the integrated conductance obtained from a 
conventional layered-earth inversion is first calculated as a 
function of depth z. The "average" or "effective" conduc- 
tivity at each depth, given by 



0"av(z) = 



S(z) 



(7) 



(Jain, 1966), is also calculated. This process is continued at 
increasing depths until z reaches the same value as that 
obtained from equation (5) at the lowest measurement fre- 
quency, / min . The value of z so obtained will be the maxi- 
mum depth of investigation z max . It should be stressed that 
even though the layered-earth linearized inversion may 
result in interpreted layers deeper than this value, then- 
existence is not supported by the data. In addition, care must 
be taken that static shifts have not corrupted the inversion 
results. This is a particular problem with the MT method and 
is discussed in more detail in the "Discussion" section. 

This relatively straightforward approach to calculating the 
depth of investigation can be tested against the rigorous but 
computationally formidable method described by Parker 
(1982). Using quadratic programming techniques, Parker 
shows that there is a "critical depth" below which the model 
conductivity can be chosen arbitrarily without degrading the 
goodness of fit of the data to the model. His critical depth is 
thus, in effect, the practical depth of investigation. 

I tested a data set from a long-period MT experiment 
described by Parker and Whaler (1981) which had been 
inverted to a number of different model solutions including 
simple layered models, smooth models, and delta function 
D + models. Although the models were very different, their 
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integrated conductance-versus-depth functions were reason- 
ably alike. One inversion result from Parker and Whaler's 
paper is shown in Figure 3a. The procedure used to calculate 
the maximum depth of investigation z max follows that out- 
lined above: first, the integrated conductance S(z) and the 
average conductivity <x av (z) are calculated from equations (6) 
and (7); these are plotted in Figures 3a and 3b. Values of the 
apparent maximum depth z msoLa (z) are then calculated for 
successive depths by substituting / = / min and <r 1 = a av (z) 
into equation (5), i.e., z maXfl (z) = 750 [a av (z)/ min ]~ 1/2 . 

The value of z for which z = z maXa is given by the 
intersection of the two lines in Figure 3c, and in this example 
equals 620 m. The certainty with which the depth of inves- 
tigation can be estimated is determined by the angle sub- 
tended by these two lines. As can be seen from Figure 3c, 
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Fig. 3 . Procedure for calculating the depth of investigation 
for a layered-earth model, (a) The model and its associated 
average conductivity cr av (z) calculated from the integrated 
conductance S(z) shown in (b). (c) shows the apparent 
maximum depth of investigation z maXo (z) calculated from 
equation (5) at the lowest measurement frequency, using 
CT i = ^avU)- The maximum depth of investigation is the value 
for which z = z max - z max . (The model is adapted from 
Parker and Whaler, 1981.) 



z maXa (z) decreases rapidly when conductive layers are e ^ 
countered and increases slowly in resistive layers. 

The depth of investigation calculated with this procedure 
is within 10 percent of Parker and Whaler's estimate 
Equally good estimates were obtained for other models ^ 
Parker and Whaler's paper, including smooth and delta- 
function models. In contrast to Parker an4- Whaler's tech, 
nique, however, the approximate scheme suggested above j$ 
very simple to implement and easy to calculate. 

Transient electromagnetic sounding 

A similar analysis can be applied to TEM data. In this 
case, the upper limit of depth of investigation is controlled 
by the earliest sample time of measurement (50 |xs in many 
systems, but as early as 6 |xs in some). Unlike the MT case, 
however, the lower limit is determined by the time at which 
the signal decays into noise. Two commonly used configu- 
rations are considered: the in-loop or central induction 
configuration in which a vertical-axis receiver is located at 
the center of a transmitter loop, and the wire-loop config- 
uration which employs a grounded wire source. 



Time of departure. — A suite of TEM magnetic-field and 
time-derivative (voltage) response curves for the in-loop 
configuration is shown in Figures 4 and 5. These data are for 
a step turnoff of current into a transmitter loop of radius a 
over a two-layer earth with upper layer thickness d. The 
vertical axis in these graphs is the measured response, either 
magnetic field h z = BJ\iq or induced voltage dBjdt\ this 
provides a more realistic estimate of resolution than does 
apparent conductivity, which is a nonlinear function. The | 
horizontal axis is presented in two forms: the bottom axis is jj 
normalized time x, and the top axis is the thickness of the | 
upper layer expressed in diffusion depths, defined by equa-|. 
tion (4). t is a dimensionless quantity which reduces to time | 
for a fixed geologic model. i 

Similar responses to these are obtained for the wire-loop | 
configuration, with a grounded-wire source of length dlt 
separated from the vertical-axis receiver by a distance a, I 
located on an axis perpendicular to the transmitter. The| 
response for this geometry can be obtained by multiplying | 
the in-loop response by dill-no. to adjust for source moment, | 
noting that a is now the source-receiver separation. The| 
ordinate becomes VlTrvcflldtm for voltage response (Figure | 
4), and h z 27ra 2 /Idi for magnetic-field response (Figure 5). | 

Consider first the two sets of curves for voltage response | 
shown in Figure 4. Values of aid = 25 represent "far-zone"! 
soundings for which the loop size or source-receiver sepa-| 
ration is much larger than the diffusion depth, and values of | 
aid = I are for 4 'near-zone' ' soundings for which the loop| 
size or source-receiver separation is less than the diffusion | 
depth. As noted by Raiche and Spies (1981), Kauahikauaf 
(1981), and Spies and Eggers (1986), the time at which the| 
second layer is first detected (termed here the time off* 
departure) is reasonably constant for all models. The depth | 
of penetration depends primarily on the sample time and noth 
the loop size or source-receiver separation. Changes $i 
conductivity at depths greater than the "depth of penetra-'ri 
tion" impart a negligible change to the response at times * 

$ 
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H'er than the time of departure. In most cases, the lower 
layer is evident at a normalized time 



It 



o> 0 $ 



- 1. 



(8) 



Another way of expressing this is that the lower layer can be 
detected beneath an upper layer with thickness equal to the 
diffusion depth 5 XD as defined by equation (4). 

The departure of the curves from the response of the top 
layer is more readily observed for some models than others. 
For example, the separation between the curves, and hence 
the detectability, increases with conductivity contrast. As 
with the MT case, the value chosen for d/h or t depends on 
the accuracy of measurement. For a 20 percent difference 
between curves with a conductivity contrast of 10, the 
departure point ranges from t = 0.4 to 4.0 for voltage 
measurements, and 0.7 to 6.0 for magnetic-field measure- 
ments. At t = 1 , the difference between the curves ranges 
from 1 percent to 80 percent, depending on the model. 

For the magnetic-field response shown in Figure 5, the 
separation between the far-zone (aid = 25) curves is gradual; 
and it is difficult to choose an appropriate value for the time 
of departure because the magnetic-field response at early 
times, given by 



6t 



ot|ul 0 a' 



(9) 



decays linearly from the primary field value at t - 0. It is 
only when intermediate times are reached (t = 10 to 20) that 
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the decay rate increases to more easily observed values. 
Although the separation between the curves at earlier times 
would be . much greater if the data were displayed as apparent 
conductivity, it would not be possible in practice to measure 
with sufficient accuracy to make these values meaningful. In 
general, for near-zone and intermediate-zone sounding, a 
conductive layer is more easily detected at a given time than 
is a resistive layer. For far-zone sounding, conductive or 
resistive layers have similar sensitivities for voltage (dBldi) 
measurements. 

Sensitivity functions for the in-loop TEM response for 
both magnetic-field and voltage measurements are plotted in 
Figure-6. These were computed numerically by perturbing 
the conductivity of a thin layer in an otherwise uniform 
half-space. The thickness of the layer was set equal to 
one-tenth of its depth and its conductivity was twice that of 
the background. The function y(z) is the percentage change 
in measured response resulting from the conductivity per- 
turbation. For the case of a thin resistive layer (having a 
conductivity contrast the reciprocal of that used for the 
conductive layer), the y(z) functions retain the same shape; 
but the amplitudes are reduced to one-half those shown and 
the signs are reversed. For near-zone sounding (aid < 1), the 
maximum contribution originates from a depth of 0.45 diffu- 
sion depths for magnetic-field measurements and 0.35 diffu- 
sion depths for voltage measurements. As expected from 
inspection of Figure 4, the undershoot for the voltage 
response increases in magnitude as the aid ratio increases. 
Approximately 10 percent of the response for near-zone 
sounding is obtained at values of dl§ > 1, i.e., when the 
buried layer is overlain by one diffusion depth of overbur- 
den. This is true both for voltage and magnetic-field mea- 
surements, except for the magnetic-field far-zone case dis- 
cussed earlier. In the following discussion, a value of t = 1 
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Fig. 5. TEM magnetic field response curves for a two-layer 
earth. The second layer is evident at dlh = 1 for the near 
zone (aid = i), but is difficult to discern in the far zone 
(aid = 25) until much later in time. 
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is taken to be a reasonable value for the time of departure for 
all but far-zone magnetic-field measurements. As is shown 
below, this value can be modified to suit the preference of 
the reader without altering the main conclusions of the 
paper. 

It is instructive to compare these curves with the corre- 
sponding sensitivity functions for magnetoteilurics, which 
are displayed as dashed lines. The MT sensitivity functions 
are based on the measured \EiH\ values, not on apparent 
conductivity, in order to facilitate comparison with the TEM 
curves. The curve plotted on the top graph is the time- 
domain MT response (Kunetz, 1972), and the curve on the 
lower graph is the conventional frequency-domain response. 
Twenty-five percent of the time-domain MT response origi- 



nates from d/b > 1 and 3 percent from dlh > 1.5. i n 



Step response 
(magnetic field) 



Y{z),% 




Fig. 6. Sensitivity functions for TEM, based on the re- 
sponse of a thin conductive layer at depth d, calculated for a 
magnetic-field receiver (top) and a voltage receiver (bottom). 
Corresponding MT sensitivity functions for \E!H I measure- 
ments (dashed) are superimposed on the graphs. The basic 
features are similar to the TEM curves, but the amplitudes 
are lower. 



frequency domain, the respective numbers are 26 percejj 
and 18 percent. Comparing the curves in Figure 6, j t . 
evident that TEM near-zone sounding {aid < 1) has ^ 
maximum sensitivity. 



Depth of investigation. — The practical limitation on 
depth of investigation for TEM systems is determined by ^ 
time t' at which the signal decays to the noise level, 
source moment, and earth resistivity. The time at which^ 
buried layer can be first detected [/ in equation (8)] j s 
function of the depth to the layer and the average resistivity 
of- the overlying section. The two times t and t' can b* 
equated to yield the thickness of overburden for which the 
TEM signal equals the noise level for any combination 
source moment and earth resistivity. The thickness of over- 
burden for which the signal equals the noise level is defined 
to be the maximum depth of investigation. 

Noise sources not geologic in origin include geomagnetic 
micropulsations, spherics, cultural noise (power lines), wind 
noise, and instrumental or amplifier noise. The recorded 
noise is drastically reduced by time averaging (grouping thi 
samples into time windows) and ensemble averaging 
(stacking repeated transients) and can be further reduced 
noise prediction techniques (Spies, 1988). In most TEM 
systems the window width is a fraction a of the sample time 
/. Typically, field measurements are averaged for 15 minutes 
or so at each station in an attempt to achieve a reasonable 
noise level. The total number of samples averaged in the 
time window immediately preceding the next current pulse is 
equal to T/2aAt, where T is the total data acquisition time! 
and At is the digitizing interval. Thus the total number off 
samples does not depend on the sample time t or the ground§ 
conductivity. In more conductive areas, the measured tran-| 
sient is longer and fewer transients are recorded in a given| 
acquisition time, but the window width increases in suchai 
way that the total number of averaged samples remains! 
constant. For statistically independent white band-limitedp 
noise, the resultant noise level after averaging is independent! 
of sample time, and the effective bandwidth BW is equal to| 
UaT. Although in practice the observed noise spectrum isg 
not white and the samples are not statistically independent,| 
there is a large enough variability to make these assumptions^ 
reasonable for the analysis which follows, especially in the! 
frequency range 0.1 Hz to 1.0 kHz applicable to deep TEM| 
soundings. It is interesting to note that remote reference! 
noise-cancellation schemes effectively whiten the resultanl| 
spectrum (Nichols et al., 1988). | 
The spectral shape and amplitude observed in practice can| 
be highly variable and depend on such factors as geographic| 
location and season (spherics are highest in the summer srif 
at mid-latitudes), time of day, and weather conditions. J| 
assume that estimates of the noise levels obtained afteif 
windowing, stacking, and other noise-reduction schemes are| 
available from field measurements or published data for the 
survey area of interest. Typically, values obtained after 15 f 
minutes of stacking range from 10~ 3 nT in winter to 2 x 10 *fr 
nT in summer for magnetic-field measurements (denoted b?^ 
t\ b ), and from 0.1 to 0.5 nV/m 2 in winter and 2 to 10 nV/m*j; 
in summer for voltage (dBIdt) measurements (denoted by : > 
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Near-zone sounding. — Near-zone, or short-offset, sound- 
refers to cases in which the loop size or source-receiver 
ing tion i s less than the depth of investigation; the induc- 
SC n number 2//a> 0 a 1 is greater than unity. The response of 



tion nuin^^ 1 — r-v ■ - - - -* - -- 

. buried layer is observed on the descending branch of the 
half-space curve (the so-called "late-time" , region), as 
! Figures 4 and 5 {aid = 1). The derivations which 



shown i» *e * 

follow are for the in-loop configuration, but can also be used 
for the wire-loop configuration in slightly modified form. 

(a) Magnetic field response. — The signal amplitude at times 
when the second layer is detected in near-zone sounding is 
given by the late-time asymptotic expression for a homoge- 
neous half-space with conductivity <x for the upper layer, 



4- 



I<r m ^ 3 0 n a 2 
30V^ 3/2 



(10) 



(Spies and Eggers, 1986). Equating the time term in equation 
(10) with the time of departure [equation (8)], it can be shown 
that 



^-0.26 



IA 

~h 7 



1/3 



(ID 



where A is the loop area and / is the transmitted current. 
Setting B z = u, 0 h z equal to the system noise level v\ B , 

1/3 



</~2.8 x 10" 3 



IA 



For a typical noise level of 10 3 nT, this becomes 



d « 28 (IA) 



1/3 



(12) 



(13) 



An example of the range of depths' investigated with mag- 
netic-field near-zone sounding is shown in Figure 7. The 
maximum depth is given by equation (12). The minimum 
depth is determined by the earliest sample time and is found 
by solving for d in equation (8). Above this depth only the 
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W i t k^ e P th of investigation (shown shaded) attainable 
noise i , ?y st< i m measuring magnetic field, assuming a 
maxim L 10 nT and lo °P parameters shown. The 
conduct™ pth fe Iven by ec * uation (12)] does not depend on 
cuvity. The minimum depth is given by equation (8). 



integrated conductance for the overlying section can be 
determined. 

(b) Voltage response. — Similarly, the late-time voltage 
(dB/dt) response for the in-loop configuration 



30\/rrf 



5(2 



and equation (8) can be combined to obtain 



d « 0.55 



I A 



1/5 



(14) 



(15) 



As before, we can equate the response V with the noise level 
t) v to obtain 



d = 0.55 



IA 



OT1 v 



1/5 



For a typical noise level t\ v of 0.5 nV/m 2 , 

f IA \ 1/5 

d~40 



(16) 



(17) 



These results are displayed graphically in Figure 8, using the 
same source parameters as Figure 7. 

The most surprising result of these derivations is that, 
unlike the result obtained for voltage measurements [equa- 
tion (15)], the depth of investigation obtained with magnetic- 
field measurements does not depend on the conductivity of 
the ground. This result appears at first to be counterintuitive. 
However, it can be seen from equation (10) that the depen- 
dence of magnetic-field response on time is the reciprocal of 
its dependence on conductivity, and thus, for any fixed value 
of h z9 an increase in conductivity is matched directly by a 
corresponding increase in sample time. Therefore, as long as 
the noise level is independent of time, the depth of investi- 



1 L_" 1 1 



Too shallow for 
TEM at 50 us 




\Q 2 

RESISTIVITY (Q m) 

Fig. 8. Depth of investigation (shown shaded) attainable 
with a TEM system measuring voltage, assuming a noise 
level of 0.5 nV/m 2 . The maximum depth [given by equation 
(16)] is proportional to the average resistivity of the overly- 
ing section, raised to the 1/5 power. The minimum depth is 
given by equation (8). 
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gation is independent of conductivity. For the voltage re- 
sponse [equation (14)], time and conductivity have different 
power-law dependencies. 

The relationships given for depth of investigation in equa- 
tions (12) and (16) exhibit different dependencies on source 
moment for magnetic-field and voltage response systems (1/3 
and 1/5, respectively). Doubling the source moment in- 
creases the depth of exploration by 26 percent for magnetic- 
field measurements but only 15 percent for voltage measure- 
ments. In order to double the depth of exploration, it is 
necessary to increase the transmitter moment by factors of 8 
and 32 for magnetic-field and voltage measurements, respec- 
tively. These formulas can also be used for the wire-loop 
configuration, by replacing the loop area A by adtll, where 
a is now the source-receiver separation and di is the dipoie 
length. 

Voltage measurements are superior for penetrating resis- 
tive areas, but magnetic-field measurements are superior for 
conductive areas. For the in-loop configuration, the earth 
conductivity a for which both modes have identical depths 
of penetration can be derived by combining equations (12) 
and (16) to obtain 



_3.1 x lO 11 ^ 3 

which for typical noise levels given earlier reduces to 

6.2 



(/A) 



(18) 



(19) 



For in-loop sounding with a 300 m loop and 25 A of current, 
magnetic-field measurements are superior for resistivities of 
less than 2800 ft ■ m. For a smaller 50 m loop with 10 A of 
current, this value is lowered to 140 ft ■ m. Thus for the 
noise level estimates used here, magnetic-field measure- 
ments are superior for all but the most resistive areas. 

Far-zone sounding. — The term "far zone" applies to 
sounding where the source-receiver separation or loop size 
is much greater than the depth of investigation: the induction 
number Itla^a 1 is less than about 0.1. 

The relevant expressions for voltage (dBJdt) measure- 
ments are. the early-time response 



31 



{jar 



for the in-loop configuration and 

3Idt 



Intra 



(20a) 



(20b) 



for the wire-loop configuration, assuming a receiver moment 
of 1 . The voltage response for far-zone sounding is inversely 
proportional to the cube of the loop size or the fourth power 
of the source-receiver separation and directly proportional 
to resistivity. Since these expressions are not functions of 
time, they cannot be combined with equation (8) in the same 
manner as for near-zone sounding. However, the greatest 
depth of investigation attainable with far-zone sounding is 
achieved at the latest time at which the far-zone approxima- 
tion is valid, i.e., 



It 



2-0.1. 



(2| 



Equating this with the time of departure [equation (8)], it Cat 
be shown that the loop size, or transmitter-receiver sepa a 
tion a, must be greater than about three times . the depth $ 
investigation. Substituting a = 3d in equations (20a) ^ 
(20b), the maximum depths of investigation are 



<i=-0.48 



cnq v 



1/3 



for in-loop sounding and 



/idty 4 

d - 0.28 



(22a, 



(22bi 



for wire-loop sounding. For values of a > d, the depth 4 
investigation varies inversely as a. 

The wire-loop configuration is often used for deep sound- 
ing (Keller et aL, 1984) because large source-receiver sep a . 
rations can be easily deployed and the signals have a 
relatively small dynamic range. Exploration to a depth of 
4 km can be achieved using a source-receiver separation 
12 km, a source length of 4 km, and a transmitter current of 
50 A for a 100 ft • m earth and 500 A for a 10 ft • m earth, 

It is not feasible to perform a magnetic-field sounding in 
the far zone because, as mentioned earlier, it is not possible 
to measure the slowly decaying field [given by equation (9)] 
with sufficient accuracy using present technology. In prac- 
tice the separation between the magnetic-field response of 
different geoelectric sections cannot be measured until inter- 
mediate times are reached. 
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Intermediate-zone sounding. — In practice a TEM sounding| 
may involve aspects of both near-zone and far-zone sound- 
ings if measurements are made over a wide enough time 
range. A suite of in-loop magnetic-field decay curves for| 
different loop sizes is shown in Figure 9. For a fixed sample | 
time, there is an optimum loop size for which the maximum jj 
signal strength is achieved; alternatively, for a particular| 
noise level and transmitter current, there is an optimum loop j 
size for measuring to late time. This occurs in the area of| 
transition between near-zone and far-zone sounding (the| 
intermediate zone) and is delineated by the straight line on| 
the figure. For example, if the noise level or the hJJ ordinate| 
is set at 0.033, then the optimum loop radius a = 10 m, since| 
this enables measurements as late in time as 2f/cru. =12, } 

The depth of investigation to be used as a constraint when| 
minimizing the transmitted current is derived as follows: g 
From Figure 9, the maximum signal level attainable at any J 
particular sample time is \, 



h 7 a 



= 0.33, 



which is measured at a normalized time 



It 

CT|X 



= 0.12 a 2 . 



(23a)? 



(23b); 
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(23*! 



(23| 



jatest time that measurements can be made is the time at 
W h the sienal h. in equation (23a) decays to the noise level 
W 'combining equations (23a) and (23b), this time is 



? = 6.5x 10 



(24) 



j in equation (24) with the time of departure 
[equation (8)], the maximum depth of investigation achieved 



Equating 



d~ 1.4 x 10" 7 — . 



(25) 



Combining equations (23a) and (25) and setting h z equal to 
the noise level ti^/jiq, the loop radius for which the maximum 
depth of penetration is achieved is 

a ** 2.9d, (26a) 

or, for a square loop of side L, 

L = a\^r = 5d, (26b) 

Thus, for a given transmitter current and noise level the 
maximum depth of exploration, given by equation (25), is 
achieved by using the loop size given in equation (26). For 
example, for a noise level of 0.01 nT, a depth of exploration 
of 100 m is achieved with a 500 m loop size with only 7 mA 
of current. These calculations show that if transmitter cur- 
rent is the limiting factor, then the maximum depth of 
investigation is achieved with very large loops sizes. 

However, transmitter current is seldom constrained to the 
low levels shown here. If the current is increased while 
keeping the loop size fixed, then equation (25) can no longer 
be used. The depth of investigation must be calculated using 
the complete expression for the half-space response, as 
illustrated in Figure 9, and equation (8). The depth of 
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investigation is increased by factors of 3.6, 2.3, and 2.15 for 
successive factors of 10 increase in the transmitted current. 
The result of using higher currents is to move down the 
descending branch of the half-space curve to the region 
where equation (12) holds, effectively changing the sounding 
to "near zone. " 

For voltage {dBldt) measurements, the analogous expres- 
sion to equation (24) for in-Io'op sounding is 



r-0.16a 1/3 )xo ( ~ 



2/3 



(27) 



Equating, this with the time of departure [equation (8)], we 
obtain 



4-0.55 



crn. 



1/3 



with 



a = Ud. 

For a square loop of side L, this becomes 



(28) 



(29a) 



L = a\/^ « 3d. (29b) 

As for the near-zone analysis, the optimum parameters for 
voltage response depend on the resistivity. For exploration 
to a depth of 200 m, the optimum loop size is 600 m. This can 
be obtained with a current level of 0.2 mA for a 100 fl - m 
earth and 24 mA for a 1 fl ■ m earth, assuming a noise level 
of 0.5 nV/m 2 . By comparing expressions (25) to (29), it can 
be seen that voltage measurements are preferable to mag- 
netic-field measurements for shallow intermediate-zone 
soundings, since a greater depth of investigation can be 
achieved with the same loop size. 

For wire-loop soundings, the depth of investigation ob- 
tained under the condition of minimum source moment is 



(idty* 

f-0.37 , 



0.37 | I , (30) 

which is obtained with a transmitter-receiver separation 



(31) 



For example, for exploration to a depth of 200 m, the 
optimum source-receiver separation is 260 m. For a source 
length of 50 m the transmitter current required is 8.4 mA for 
a 100 fl * m earth and 800 mA for a 1 fl - m earth. 

By comparing equations (22b) and (30), it can be seen that 
the depth of investigation obtained with far- zone sounding is 
inferior to that attainable with intermediate-zone sounding. 

Unlike the far-zone case, magnetic- field measurements 
can be made for grounded-source systems in the intermedi- 
ate zone. The relevant expressions are 



and 



a - \Jd. 



1/2 



(32) 



(33) 



The depth of investigation is again independent of conduc- 
tivity, and the dependence on transmitter moment has in- 
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creased from 1/4 to 1/2. For the same example given earlier 
(dipole length 50 m), the depth of investigation of 500 m is 
achieved by increasing the separation to 850 m and the 
current to 500 mA. For the same source moment, the depths 
of investigation attainable with magnetic-field measurements 
are much greater than for voltage measurements. 

With intermediate-zone soundings the depth of investiga- 
tion is maximized for any given transmitter current (in-loop) 
or source moment (wire-loop). In practice, intermediate- 
zone soundings are not often used for the in-loop configura- 
tion, since it is more practical to use smaller loop sizes and 
higher current. In these cases the depth of investigation is 
greater than the loop size, and the sounding reverts to being 
"near zone." 

COMPARISON OF DIFFERENT METHODS 

Although the time of departure or thickness of the over- 
lying section expressed in diffusion depths is reasonably 
constant for all systems, there are large differences in the 
maximum depths of investigation achieved in practice when 
signal and noise levels are considered. Table 1 presents a 
convenient summary of the expressions derived in this paper 
for a wide range of source-receiver configurations. 

MT is classified as "far zone" because of the great 
distance of the source. Expressions for CSAMT can be 
derived in a similar manner to that described for TEM. The 
depth of investigation of the CSAMT method differs from 



MT, since signal and noise levels need to be taken im 
account. 

The TEM expressions vary greatly for "near" and "f ar , 
zones and differ for voltage (dB/dt) and magnetic -field me* 
surements. The optimum separation and configuration de 
pend on the required depth of exploration and the averag, 
resistivity of the overlying section. 

The expressions given in Table 1 give the depth 0 ( 
investigation based on detection of a buried feature at depth 
It is usually necessary to use a larger source moment than 
that obtained from the table so that measurements can be 
made at a significantly later time or lower frequency t 0 
obtain an unambiguous signature for adequate resolution. 
The amount of additional data required in practice depends 
on such factors as system accuracy, geologic noise, the 
dimensions of the target, and its conductivity contrast with 
the overlying strata. From field experience, however, an 
additional one-third decade (i.e., a factor of two) of time or 
frequency is usually adequate. 

Applying the one-third decade guideline to the near-zone 
and intermediate- zone configurations in Table 1 reduces the 
effective depth of investigation to 71 percent of the values 
shown. The factor shown in brackets is the increase j n 
source moment required to extend measurements to twice 
the sample time. A different criterion is applied to near- zone 
soundings because, as can be seen from Figure 4, the signal 
level for a homogeneous earth is constant with time and! 
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Table 1. Depth of investigation of various EM sounding methods. 



Near zone 
It 

- 2 >5 



MT 



TEM, in-loop, voltage 



TEM, wire-loop, voltage 



TEM, in-loop, magnetic field 



TEM, wire-loop, magnetic field 



N/A 



0.55 



0.48 



I A 



(XT), 



ladt 



1/5 



[5.7] 



1/5 



[5.7] 



2.8 x 10 " 3 I — I [2.8] 



2.4 x 10" 3 



(ladt 



1/3 



[2.8] 



Intermediate zone 
(minimize current) 



N/A 

/ \ 1/3 

0.55 | [2.8] 

\<rr) v / 
when a= \Jd 

(Idt\ m 
0.37 [4] 

when a = 1 3d 

1.4 x 10 ~ 7 — [1.4] 

when a = 23d 

10 u 1/2 [2] 

when a - \.ld 



Far zone | 

It \ I 
2 ^ 0.1; a>3d \i 



750 



1/2 



1/3 



0.48 ( 1 [2] 

for a = 3d 

0.28 [2] 

for a = 3d 

can't measure 



can't measure 



ti v = voltage noise level after stacking (typically 0.5 nV/m 2 ) 
t) b - magnetic-field noise level after stacking (typically 10 " 3 nT) 

a - loop radius (in-loop) or transmitter-receiver separation (wire-loop) 

A = loop area (in-loop) 

I = transmitter current 
dt ~ grounded dipole length (wire-loop) 

22cnuK Zr dete u i0n ° f 2 buried feattirc at depth. It is usually advisable to use a larger source 

moment (multiply Id€ or 7 A by the factor shown m brackets) to facilitate interpretation (see text). 
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s w hen a conductive lower layer is detected. The 
omfliendation in this case is to double the source mo- 

m The range of depths covered in a typical MT and TEM 
1 e) s ounding is shown in Figure 10. In very conductive 
ons the TEM sounding covers several decades of depth; 
h t as the resistivity increases, the depth range narrows. As 
the loop size is reduced and the resistivity of the ground is 
increased, the depth range finally decreases to zero; at this 
ooint the signal has decayed to the noise level before the first 
Sample time. For this example there is overlap in the depth 
of investigation obtained with the two methods for resistiv- 
ities less than 260 fi • m. Adequate overlap is essential for 
successful implementation of MT static-shift corrections 
using TEM. 

An in-loop TEM (voltage) sounding has the same depth of 
investigation as MT at a frequency / if its source moment I A 
is set equal to 



I A = 



2.4 x 10 6 



V2 



f 



5/2 



(34) 



assuming a TEM noise level y\ v of 0.5 nV/m 2 . For/ = 100 Hz 
and / = 10 A, the depth of investigation in a 1 O • m earth is 
75 m (using a TEM loop size of 1.5 m), whereas a depth of 
2400 m is achieved for a 1000 ft ■ m earth (280 m loop size). 
Similar calculations can be performed for other configura- 
tions. It should be mentioned here that it is inadvisable to 
use very small TEM loops for making MT static-shift cor- 
rections because of the possibility of recording anomalous 
superparamagnetic effects (Spies and Frischknecht, 1989). 

EXAMPLE 

It is useful to work through a case history of survey design 
for a typical field problem to illustrate the use of the 




RESISTIVITY (Qm) 



with viv Comparison of regions of investigation obtained 
" n and in-loop TEM (voltage response). For this 
^ample the upper MT frequency is 100 Hz, the TEM loop 

invesr ^° m ' and the current is 10 A * The de P th s of 
than ?firfn° n of the two metn ods overlap for resistivities less 
nece* ' m * ^° ensure overlap at higher resistivities, it is 
emnil Sary i e ^ tner t0 me asure at higher MT frequencies or to 
l P l oy a larger TEM transmitter moment. 



expressions given in this paper. Consider, for example, that 
the geologic objective is to map the thickness of volcanics 
with an average resistivity of 100 H • m overlying a 5 H * m 
sedimentary section. The thickness of the volcanics is ex- 
pected to vary from 0 to 2000 m. A contractor can supply a 
TEM system capable of in-loop sounding with a transmitter 
loop size of either 100 m (10 A of current) or 400 m (20 A). 
The earliest time channel is 50 jxs. Moderate noise levels 
(0.5 nV/m 2 ) are predicted. 

Equation (8) sets the limit on the thinnest volcanic cover 
that can be mapped which, based on the earliest sample time 
of 50 |xs and a resistivity of 100 ft • m, is calculated to be 
90 m (Figure 8 is applicable to this example). Depths 
shallower than 90 m cannot be resolved with this system. 

The maximum depth of investigation is calculated from 
Table 1 using the near-zone formula, since the depth of 
investigation is expected to be greater than the loop size. 
The maximum depth at which the sedimentary basement can 
be detected is calculated to be 1000 m for the 100 m loop and 
2000 m for the 400 m loop. For a confident interpretation, 
though, these depth estimates should be reduced to 71 
percent of these values, i.e., 700 m and 1400 m. The latest 
sample time at which the signal is above the noise level 
calculated from equation (8) is 13 ms for the 100 m loop and 
53 ms for the 400 m loop. These times are used to choose the 
appropriate pulse repetition rate. To resolve an interface at a 
depth of 2000 m, the source moment (area-current product) 
must be increased by a factor of 5.7, for example, by 
employing a loop size of 800 m with 30 A of current (a source 
moment of 1.8 x 10 7 Am 2 ). 

If the system were capable of transmitting into a grounded 
source, a wire-loop sounding could be made. For far-zone or 
long-offset sounding the minimum source-receiver separa- 
tion is 6000 m (aid = 3). This would require a source moment 
Idi of 13 000 Am, e.g., 26 A transmitted into a 500 m 
grounded dipole (or 26 000 Am to allow for interpretation of 
the conductive lower layer). The same source-receiver sep- 
aration could be used for all depths between 90 m and 
2000 m, but unlike the near-zone case, the signal level would 
not increase for shallower depths of the interface. The 
smallest possible source moment (4200 Am) is achieved 
using a separation of 2600 m (intermediate- zone sounding 
with aid = 1.3). The source moment should be increased by 
a factor of 4 to 16 800 Am for interpretation. Larger source 
moments would be required for smaller separations (near- 
zone sounding), e.g., 7000 Am for a source-receiver separa- 
tion a = 0.45J = 900 m. This moment should be increased by 
a factor of 5.7 to 40 000 Am for interpretation. 

The same depth range (90 m.to 2000 m) could also be 
probed with an MT system having a frequency range 7 Hz to 
7 kHz (to facilitate interpretation, the 7 Hz value is half that 
calculated from Table 1). 

DISCUSSION 

The expressions for the depth of investigation derived for 
TEM in this paper are based on the detection of a buried 
half-space. It is logical to question whether the same con- 
clusions would apply to the detection of more general 2-D 
and 3-D bodies. Newman et al. (1986, 1987) and Newman 
and Hohmann (1988) published a number of 3-D model 
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results for the TEM in-loop configuration, and Gunderson et 
al. (1986) published results for a grounded source. All these 
models, whether a basement uplift or conductive prism, are 
first detected under approximately the same thickness of 
overburden as for the case of the buried half-space. Just as 
for a buried layer, the exact time at which an interpreter 
would be confident that a 3-D body is being detected 
depends on its size and conductivity contrast and is some- 
what subjective. For example, Eaton and Hohmann (1989), 
in a study of 3-D bodies in a layered earth, prefer an estimate 
of 1.24 diffusion depths. This would have the effect of 
reducing the depth of investigation given by equation (16) by 
11 percent. The exact value chosen will not alter the main 
conclusions reached in this paper; it will merely result in a 
slight change in the multiplicative constant. 

A basic assumption made in this paper is that the earth can 
be adequately represented by a 1-D model. This assumption 
is reasonable for most controlled-source EM soundings in 
quasi-layered environments. Recent papers (Ranganayaki, 
1984; Ingham, 1988) have demonstrated the utility of 1-D 
interpretation of MT data obtained using rotationally invari- 
ant MT parameters, as proposed by Berdichevsky and 
Dimitriev (1976a). A significant hazard with this approach, 
however, is the possibility that the data have been corrupted 
by static effects caused by near-surface inhomogeneities or 
topography (Berdichevsky and Dimitriev, 1976b; Park, 1985) 
or large-scale structure which distorts the regional current 
system (Ranganayaki and Madden, 1980). If, however, the 
soundings are corrected for static shift using schemes such 
as those proposed by Rokityansky (1982), Jones (1988), 
Sternberg et al. (1988), and Bahr (1988) for MT, or by 
Newman (1989) for TEM sounding, then 1-D inversions will 
be relatively undistorted, and the relations derived in this 
paper will give satisfactory results. 

Other sources and configurations. — Spies and Frischknecht 
(1989) show that the time of departure for TEM given by 
equation (8) is applicable to virtually any source-receiver 
combination. For grounded dipole-dipole configurations, how- 
ever, the response curves first need to be adjusted for static 
shifts caused by dc current flow in the ground. 

Similar analyses, relating the depth of investigation to 
signal and noise levels, can be applied to virtually any EM 
sounding technique. For example, Keller (1971) calculated 
similar relationships for a frequency-domain loop-loop sys- 
tem used for sounding in Hawaii. He concluded that the 
optimum source-receiver separation is twice the depth of 
investigation, and that the depth of investigation is propor- 
tional to the source moment raised to the 2/9 power. Of 
course the most complete analysis for any EM system is to 
calculate the EM response over the full range of field 
parameters, taking into account constraints imposed by 
signal and noise levels, and measurement accuracy. 

Frequency-domain controlled-source methods using a 
loop-loop (slingram) configuration are often used for shallow 
sounding, since the operating frequency can be made arbi- 
trarily high and the skin depth reduced accordingly. How- 
ever, at very high frequencies displacement currents can be 
significant, thereby complicating interpretation. An alterna- 
tive approach for shallow sounding is to measure the quadra- 
ture response at very low induction numbers, where the 



penetration of the fields is not limited by the skin effect, 
rather by the source-receiver separation r. It can be shov^ 
that approximately 75 percent of the response at low induc- 
tion numbers originates from depths of less than % 
(Frischknecht et al., 1989). Most portable terrain conductiv- 
ity meters are based on this principle. 

CONCLUSIONS 

The depth of investigation of electromagnetic sounding 
methods is determined by the frequency or time of measure, 
ment, the earth conductivity, and, for controlled-source 
methods, the source moment and noise levels. The most 
notable difference among various methods is their sensitivity 
to the resistivity of the overlying section. The methods 
considered here, numbered in order of increasing depen- 
dence on resistivity, are: 

(1) TEM (magnetic field) ... no dependence 

(2) TEM (wire-loop and in-loop, voltage, near zone) . . . 
depth cc p 1/5 

(3) TEM (wire-loop, voltage, far zone) . . . depth oc p 1A » 

(4) TEM (in loop, voltage, far zone) . . . depth cc p 1/3 1 

(5) MT . . . depth a p 1/2 I 

The MT method, with appropriate static-shift corrections, I 
should be used wherever very large depths of investigation 
are required. For intermediate depths there is overlap be- 
tween the MT and TEM methods, which, except for TEM 
magnetic-field measurements, depends strongly on the earth 
conductivity. 

For a layered-earth section the depth of investigation 
depends on the average conductivity (defined as the inte- | 
grated conductance divided by the total thickness) of the \ 
overlying section. The average conductivity, calculated as a I 
function of depth, can also be used. for estimating an approx- § 
imate lower depth bound of a layered-earth inversion. The | 
lower bound is the depth which is equal to 1 .5 skin depths in | 
the frequency domain, or 1 diffusion depth in the time | 
domain, calculated using the average conductivity at the | 
lowest measurement frequency or latest sample time. f 

The highest measurement frequency or earliest sample | 
time determines the minimum depth at which individual I 
conductivity variations can be resolved. At depths shallower jjj 
than approximately 1.5 skin depths or 1 diffusion depth, | 
calculated at the highest measurement frequency or earliest | 
sample time, the response is determined mainly by the | 
average conductivity of the overlying section. \ 

For applications where resolution at very shallow depths |: 
is required, such as groundwater, engineering, environmen- | 
tal, and archaeological applications, frequency-domain sys- g 
terns offer the best solution. The earliest sample time in most | 
TEM systems is about 50 ^s, which is unsuitable for shallow \ 
applications. ? 
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APPENDIX A 
CONDUCTIVITY IMAGING — A REVIEW 



Most approximate techniques for interpreting EM sound- 
ing data are based on direct transformation of observed EM 
data to inferred conductivity and depth. These techniques, 
which can be collectively referred to as "conductivity imag- 
ing' * methods, demonstrate the relationship between skin 
depth (or diffusion depth) and the effective depth of penetra- 
tion of EM fields in the earth. 

An approximate depth transformation commonly used in 
magnetotellurics is known as the Bostick transform (Bos- 
tick, 1977) or the Niblett Sayn-Wittgenstein transform (Nib- 
Iett and Sayn-Wittgenstein, 1960). The Bostick transform is 
based on the asymptotic solutions for a two-layered earth 
with an infinitely conductive or resistive basement, and the 
Niblett Sayn-Wittgenstein transform is derived directly from 
Maxwell's equations; these two transforms can be shown, 
however, to be equivalent (Jones, 1983). The effective depth 
of penetration h for these transformations estimated directly 
from the Cagniard apparent resistivity p„(w) is defined by 

/p *(<■>) 

h= J , (A-l) 

which, for a homogeneous earth, is l/V2or 71 percent of the 
skin depth given by equation (2). Schmucker (1970, 1979) 
presents a related approach in which the earth is modeled by 
a "perfect substitute conductor" at a depth h*. Weidelt 
(1972) refers to this depth as the "center of gravity of the 
in-phase induced current system." The relationship between 
Bostick's and Schmucker's transformations is discussed by 
Schmucker (1987); for a homogeneous earth the Schmucker 
depth is equal to 71 percent of the Bostick-Niblett depth, or 
0.5 skin depths. The depth given by equation (A-l) is also 



used in modified forms of the Bostick-Niblett transform 
described by Jones and Foster (1986) and Murakami (1985). I 
Estimates for depth of penetration used in other imaging I 
schemes sometimes differ from the figure of 0.71 skin depths | 
used above: Jain (1966) defines his depth of penetration as f 
0.67 skin depths, and Gamble (1983) uses a figure of 0.62 skin | 
depths. The "depth of penetration" used in this context r 
refers to the depths associated with the corresponding inter- I 
preted conductivities, and, as shown earlier, considerably | 
underestimates the actual depth of investigation defined jf 
using detection criteria. ) 

Related schemes have also been used for controlled- | 
source FEM data in recent years. Sinha (1979), Wilt and 
Stark (1982), and Mundry and Blohm (1987) describe ap- 
proximate interpretation schemes for loop-loop FEM sound- 
ing data which are essentially identical to the Bostick-Niblett 
transform. 

Another approach, most commonly applied to TEM data, 
is based on the calculation of "image depths" in the earth 
from the inferred diffusion velocity. This originates with 
Maxwell's (1892) observation that the time-domain fields 
observed above a thin conductive sheet following a step 
change in magnetic field are identical to those generated by 
an image of the source receding with a velocity inversely 
proportional to the conductance of the sheet. However, the 
image location for a thin sheet bears no relation to the true 
current distribution in the sheet. Nabighian (1979) developed 
an analogous concept for a homogeneous earth and showed 
that the observed effect of currents diffusing into the earth 
from a loop source could be approximated by a single 
current filament moving downward with a velocity 
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j nC reasin2 it s horizo ntal dimension in proportion to the 
on depth V2 r/o> 0 . This filament moves downward 
doutward from the loop in a fashion similar to the actual 
a " . cur rents, but with a higher velocity and at a steeper 
I of 47 degrees compared to the 30-degree angle for the 
Uus of the maximum induced currents. 
° These ideas were used by Macnae and Lamontagne (1987) 
who developed an approximate conductivity-depth imaging 
I orithm based on inferred image depths for measurements 
made both inside and outside a large transmitter loop. They 
used a combination of 14 discrete images centered at a 
"reference depth" h to show that the slowness dtldh of the 
center of the image distribution was very nearly proportional 
to the cumulative conductance measured from the surface 
down to the reference depth and could thus be used to 
estimate conductivity. Their penetration depth was set equal 
to 54 percent of the image depth. Nekut (1987) and Eaton 
and Hohmann (1989) simplified this approach by using a 



single image to construct the approximate earth section. 
Eaton and Hohmann found that the best depth estimates 
were 44 percent of the image depth for a loop or grounded 
source, and vertical-axis receiver and 33 percent for a 
grounded source and horizontal-axis receiver. A minor ad- 
justment to these, depth estimates was made at early times. 
Nekut found that it was necessary to modify the depth 
according to the induction number. For large loop sizes (or 
early times) the depth was set at 71 percent of the image 
depth, and for small loop sizes (or late times) it was set at 48 
percent. These adjustments are not necessary in Macnae and 
Lamontagne's scheme because the width of their image 
distribution increases with time. 

Frequency-domain solutions using complex images have 
been used by numerous authors for computing rapid approx- 
imate solutions in EM induction studies. Thomson and 
Weaver (1975) give solutions for an arbitrary magnetic 
source above a muitilayered earth, and Wait and Spies (1969) 
derive image solutions for a line current source. However, 
interpretation schemes based on image depths are not widely 
used for interpreting FEM field data. 



APPENDIX B 
FRECHET DERIVATIVES FOR A PLANE- WAVE SOURCE 
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The Frechet derivatives for a homogeneous half-space 
subject to plane- wave excitation, based on an apparent 
conductivity definition, are: 

Act,,/ («) = p 2 Re {ye " 2 ^> Aa(z) dz 
Jo 

for the frequency domain, 



Jo 



to aJ {t) = I 8zd 2 e~ 4zV ' &<r(z)dz 



for the time-domain step response, and 



to aJ (t) 



Jo 



8 z$ 2 e 
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(3 - 8 z l B z ) Act(z) dz 



for the time -domain impulse response, where 7 = V/o>o> 0 , 
and 9 = Vap, 0 /4r, respectively (Gomez-Trevino, 1987b). 
to a is the change in observed apparent conductivity result- 

from a small change Aa(z) in the conductivity function 
°iz)- The term within the integral sign preceding koiz)dz is 
the Frechet derivative relating the small perturbations Aa(z) 
and Aa„. The Frechet derivatives or sensitivity functions 
also play the role of weighting functions operating on the 
conductivity distribution itself, from which the observed 
responses can be synthesized (Gomez-Trevino, 1987a). Each 

he a a values represents a distinct weighted average of 
pjr" Tne corresponding weighting functions, denoted by 

^ are plotted in Figure B-L 
of th* °^ served a PParent conductivity is a weighted average 

he true conductivity function a(z) from z = 0 to a certain 



depth z max . Although there is theoretically a finite contribu- 
tion from all depths, the value of z max depends principally on 
the system accuracy (Weidelt, 1972; Fischer et al., 1981). 
The step response gives well-behaved localized averages at 
depth, with the maximum contribution, and hence sensitiv 



T.D. (impulse) 




Fig. B-l. Apparent conductivity sensitivity functions for 
plane-wave exitation of a homogeneous earth in both the 
frequency domain (F.D.) and time domain (T.D.). Apparent 
conductivity is a weighted average of the actual earth 
conductivity from the surface down to a limiting depth 
controlled by the system accuracy. About 5 percent of the 
response for F.D. and T.D. (impulse) response, and 1 
percent of the step response, originates from depths greater 
than 1.5 5. Here 8 is the skin depth in the frequency 
domain and the diffusion depth 8 TO in the time domain (after 
Gomez-Trevino, 1987b). 
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ity, at half the diffusion depth S TD (Figure B-l). Based on the 
area under the curves, 87 percent of the response originates 
from depths shallower than 8 TD , and 99 percent from depths 
less than 1.5 8 TD . The frequency-domain and the time- 
domain impulse responses both have side lobes which are 
manifested as oscillations in apparent conductivity curves at 
conductivity discontinuitues. For the time-domain impulse 
response, the maximum sensitivity is at 0.37 8 TD . There is a 
positive contribution from depths less than 0.86 8^, and a 



negative contribution from greater depths; 78 percent of the 
total response originates from depths less than 8-^, and 95 
percent from less than 1 .5 S TD . For the frequency-domain 
case, 95 percent of the response originates from depths 
shallower than 1.5 8 TD . Since the shape of the sensitivity 
function depends on the definition used for apparent conduc- 
tivity (e.g., Spies and Eggers, 1986), caution should be 
exercised in attaching too much significance to the absolute 
percentage values given here. 
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